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Abstract. The quantum theory of the cold atom micromaser including the effects of gravity is established 
in the general case where the cavity mode and the atomic transition frequencies are detuned. We show that 
atoms which classically would not reach the interaction region are able to emit a photon inside the cavity. 
The system turns out to be extremely sensitive to the detuning and in particular to its sign. A method 
to solve the equations of motion for non resonant atom-field interaction and arbitrary cavity modes is 
presented. 
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The coupling of atomic motion to light is a topic of ma- 
jor interest in quantum optics. The mechanical effects of 
(laser) light on atoms can be exploited to cool, trap and 
handle the atoms with a great accuracy [1] (see for ex- 
ample [2] for manipulation of individual atoms with op- 
tical tweezers). The achieved control on the atomic mo- 
tion gives rise to a host of applications, like the develop- 
ment of new optical frequency standards based on single 
ions [3J, the realization of an atom laser 0], the possi- 
bility of implementing electronic components with atoms 
(instead of electrons) [SBl] , or quantum information pro- 
cessing with cold atoms and trapped ions [8]. In the last 
example, the interplay between external (motional) and 
internal degrees of freedom plays a prominent role. This 
is also the case in the usual micromaser (see Fig. [T]) when 
cold atoms rather than thermal ones are sent through the 
cavity. Indeed, in this regime a new type of induced emis- 
sion has been shown to occur because the wave behavior 
of the atoms becomes important [9]. The resulting pro- 
cess called microwave amplification via z-motion-induced 
emission of radiation (mazer) [9] has opened up a new 
chapter of micromaser physics pTT |[TTVT^1Jl8iri^ fl% |riWT^ 
fT8lfT9j . In all these papers, gravity effects on the quan- 
tized atomic motion are not considered. However, it has 
often been argued that the achievement of the cold atom 
micromaser would be a formidable experimental task as 
the atoms in this regime move so slowly that they start 
exactly to be extremely sensitive to earth gravity and are 
expected to fall down before entering or leaving the cavity 
region, breaking the unidirectionality of the atomic tra- 
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Fig. 1. Scheme of the horizontal micromaser. 



jectories [20 . A possible way to deal with gravity for the 
mazer action while keeping the atomic trajectories one- 
dimensional is to consider a vertical geometry where the 
atoms are sent vertically in the direction of the cavity, 
similarly to cavity QED experiments reported in [21,22, 
23,24,25,26,27,28]. Recently, we established the quantum 
theory of such a vertical mazer, taking into account grav- 
ity effects on the vertical quantized atomic motion |29j . 
The theory was written for two-level atoms in the reso- 
nant case where the cavity mode frequency u> is equal to 
the atomic transition frequency luq. Here, we remove this 
restriction and investigate the effects of a detuning on the 
quantum evolution of the combined atom-cavity system in 
a vertical configuration where gravity action is taken into 
account. 

The paper is organized as follows. In Sec. [2j the Hamil- 
tonian modeling our system is presented. The properties 
of the induced emission probability of a photon inside the 
cavity are then presented in Sec. [3j The connection with 
the classical regime is discussed. A brief summary of our 
results is finally given in Sec. [4j In Appendix A, a tech- 
nical derivation of a specific formula used in our paper is 
derived. The method we have developed to solve numer- 
ically the equations of motion of the mazer in the non 
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Fig. 2. Scheme of the verticai mazer. 



resonant case for any mode function is given in Appendix 
B. 

2 Model 

2.1 The Hamiltonian 

We consider a two-level atom moving in the gravity field 
along the vertical z direction on the way to a cavity which 
we define to be located in the range < z < L (see 
Fig. [J). No transverse motion is considered. The atom 
is coupled nonresonantly to a single mode of the quan- 
tized field present in the cavity. The atomic center-of-mass 
motion is described quantum mechanically and the usual 
rotating-wave approximation is made. The Hamiltonian 
thus reads 



V 

2m 



mQz 



hgu{z){a) a + acr'), 
(1) 

where p is the atomic center-of-mass momentum along 
the z axis, m is the atomic mass, Q is the acceleration 
due to gravity, luq is the atomic transition frequency, lu 
is the cavity field mode frequency, a — \b)(a\ (\a) and \b) 
are, respectively, the upper and lower levels of the two- 
level atom), a and a' are, respectively, the annihilation 
and creation operators of the cavity radiation field, g is 
the atom-field coupling strength, and u(z) the cavity field 
mode function. We denote in the following the detuning 



u> — ujq by 5, the cavity field eigenstates by |n), the global 
state of the atom-field system at any time t by \ip(t)), and 
the classical height attained by an atom of energy E in 
the gravity field by 



E 
mQ 



(2) 



2.2 The wave functions 



We introduce the orthonormal basis 



\r+(0)) = cos9\a,n) +sm9\b,n + l), 
|r„(«)) = -sin0|a,ra) + cos<9 \b, n + 1) 



(3) 



with 9 an arbitrary parameter. The \r^{9)) states coin- 
cide with the noncoupled states |a, n) and |6, n + 1) when 
9 = and with the dressed states when 9 = 9 n given by 



cot26>„ =— — - 
j / n 

with the Rabi frequency 

fl n = 2g\/n + 1 



(4) 



(5) 



We denote as |±,n) the dressed states \r^(9 n )). The 
Schrodinger equation reads in the z representation and in 
the basis ([3]) 



d 



2m dz 2 



+ mQz + (n + l)Hui 



cos 9HS + hg «(z)V«+l sin 29 
sin 29 



hgu(z)\/n+l cos 29 



2 9h5- 



2m dz 



hgu(z)y/n + l cos 29 



hgu(z)y/n + l sin 29 
sin 29 



mQz + (n + l)hu 



(6a) 



(6b) 



ns 



ith 



i>±{z,t) = {z,r±{e)m))- 



(7) 



We get for each n two coupled partial differential equa- 
tions. In the resonant case (<5 = 0), these equations may 
be decoupled over the entire z axis when working in the 
dressed state basis and the atom-field interaction reduces 
to an elementary scattering problem in the presence of the 
gravitational field over a potential barrier and a potential 
well defined by the cavity [29J. In the presence of a de- 
tuning, this is no longer the case : there is no basis where 
Eqs. ([6]) would separate over the entire z axis and the in- 
terpretation of the atomic interaction with the cavity as a 
scattering problem over two potentials is less obvious as it 
is also the case when gravity is not taken into account |16j . 
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In the noncoupled state basis (6 = 0), Eqs. ([6|) read 



2m dz 2 



mQz 



r n {z,t) 



(8a) 



+ Hgu(z)^n+1 ip n+1 (z,t), 



"at'- 



h2 d2 n *r 

■ ~z + myz + Fid 

2m oz z 



+ hgu(z)y/ n+I ip„(z,t), 



with 



i> b n+1 (z,t) = j(<*>+™)* ( Zjfe>n + i|^( t )). 



(8b) 

(9a) 
(9b) 



In Eqs. (J9j> , we have introduced the exponential fac- 
tor e l ( u a+nw)t m or der to define the energy scale origin 
at the | a, n) level. If we assume initially a monoenergetic 
excited atom coming upwards upon the cavity that con- 
tains n photons, the atom-field system is described outside 
the cavity by the wave function components (which cor- 
respond to the eigenstate \4>e) of energy E) 



_ e -iEt/H ( 



r n {z,t) = e- tmin V a E , n {z) 

^n+lO 2 ' t) = e r L . „.. 



(10a) 
(10b) 



with <^l : n( z ) and V E ,n+i{ z ) obeying 
h 2 d 2 



2m dz 2 
2m dz 2 



+ mQz-E\ip E (z) 



Hgu(z)Vn+l <p E , n+1 (z), (11a) 



(lib) 



iGz-E + M)ip Etn+1 (z) 



able 
with 



hgu(z)Vn+l(p% in (z). 
Introducing, for any length x, the dimensionless vari- 

x = x/e (12) 

1/3 



2m 2 g 



(13) 



Eqs. dTTJ) read 



d^ 

dz 2 



z + h E 
z + h E 



<P a E ,n{z) = hint u(z) <p b E>n+ i(z), (14a) 



& ) <PE, n +l{2) = h int u(z) <p%Jz), 



where 
and 



hint = gVn + 1, 

hg ~ H5 
mQl mQl 



(14b) 
(15) 
(16) 



The characteristic length £ (Eq. (113[1 ) gives the spatial 
scale of the oscillations of the wave function of a particle in 
the gravitational field near its turning point. For rubidium 
atoms, £ ~ 0.3 /im and, roughly, the dimensionless variable 
5 is numerically equal to the detuning 5/2n expressed in 
kHz. Similarly h E and L yield the classical height h E and 
the cavity length L, respectively, in thirds of /urn [29] . 

Outside the cavity, the mode function u(z) vanishes 
and Eqs. fli]) read 



(17a) 
(17b) 



— -z + h E ) y% n (z) =0, 



d 2 

- z + h E - 5 ) <y5£ : „+i(5) = 0. 



The solutions to Eqs. (|17[) are given by linear combina- 
tions of Airy Ai and Bi functions [30] and may be written 
in the form 



Ai(z — He) 



z > L 



\J(z-h E ) + d% n ~D(z-h E ) z<0 




n+1 Ai(z -h E + 5) z > L 
d b E n+1 D(z - h E + 6) z<0 



U(z) = Ai(z)+iBi(z) 
D(z) = Ai(z) -iBi(z) 



(18) 



(19) 



(20a) 
(20b) 



Equations (fl8|) and (fT9|) underline that the wave func- 
tions below the cavity consist of an upward and a down- 
ward wave (U and D, respectively). Indeed, the probability 
current densities associated to these two waves are respec- 
tively given by 

^ irm£ ^ 

For large negative z — h E values, these two waves be- 
have similarly to plane waves with a z dependent wave 
vector [29] . 

Consequently, the following interpretation must be 
given to solutions (fl"0|) . The excited atom coming upwards 
upon the cavity will be found propagating downwards in 
the upper state or in the lower state with amplitude d E 

and d E n+1 , respectively. However, in contrast to the res- 
onant case, the atom propagating downwards in the lower 
state | b) possesses a total external energy E — hS different 
from its initial value E. The atomic transition \a) — > \b) 
induced by the cavity is responsible for a change of the 
total external atomic energy. According to the sign of the 
detuning, the cavity will either provide energy to the atom 
(for 5 < 0) or remove energy from the atom (for 5 > 0). 
This results merely from energy conservation. When, af- 
ter leaving the cavity region, the atom is passed from the 
excited state \a) to the lower state \b), the photon number 
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has increased by one unit in the cavity and the internal 
energy of the atom-field system has varied by the quantity 
fku — hujQ = HS. This variation needs to be exactly coun- 
terbalanced by the external energy of the system, i.e., the 
total external atomic energy. 

Inside the cavity, the problem is much more complex 
since we have two coupled partial differential equations. 
In the special case of the mesa mode function \u(z) = 1 
inside the cavity, elsewhere], the problem is, however, 
greatly simplified. In the dressed state basis (0 = n ), the 
Schrodinger equations © take the following form inside 
the cavity : 



d 



h 2 d 2 
2m dz 2 



mQz + Vl d 



^M) (22) 



with 



4>t(z,t)=e^ +n ^(z,±,nm)) 



(23) 



and 



V+ = sin 2 n hS + hgVnTT sin 26 n , (24a) 
V- =HS- V+. (24b) 

Using Eq. Q, we have the well known relations [31j 

t&nOn = 



smt 



l mi 



Q'„ -8 ' 



COS n 



(25) 



— , cot0„ = 



n' n -s 

O'+S' 



with the generalized Rabi frequency 



We thus have 



V* = hg\J n + 1 tan n 



V n = —hgy/n + 1 cotO n . 



(26) 

(27a) 
(27b) 



The exponential factor e *(^o+n^)t nas b een introduced 
as well in Eq. (|23| in order to define the same energy scale 
inside and outside the cavity. The positive internal energy 
increases with positive detunings and vice versa with 
negative ones. For large positive (resp. negative) detun- 
ings, tends to the \b,n + 1) (resp. \a,n)) state energy. 

The most general solution of Eqs. (|22|) is given by 



^t{z,t) = e-^' h V ± E (z) 



(28) 



with 



V E M) = A± M(z-h E +h±)+B±Bi(z-h E + h±), (29) 



where A^ et are complex coefficients and 
V^/mQ, i.e. 



hn = Knt tan n , 

h~ = -h int COt n , 



(30a) 
(30b) 



From Eq. ([3]), we may express the wave function com- 
ponents of the atom-field state inside the cavity over the 
non-coupled state basis. We have 

ipn(z,t) = cos0 n ^+OM) - sm0 n ipn(z,t), (31a) 
ip n+1 (z,t) = sin0 n '4> n '(z,t) + cosO n ipn(z,t). (31b) 

This allows us to find the wave function components 
of the eigenstate \4>e) over the entire z axis. The relations 
(fTOll hold with 



a a En Ai(z - h E ) z > L, 

\J(z-h E ) + d% n D(z-h E ) z<0, 



(32) 



a E,n+i - h E + 6) z > L, 



d E,n+l V(z~h E 



z < 0, 



(33) 



and 



cos n \A+ Ai(z-h E + h n ~) 
— sin n (^A~ Ai(z — h E + h~ 



B+Bi(z-h E + h+) 
+ B~Bi(z-h E + h 



(34a) 



fE,n+i( z )\ c = 

sin On (A+ Ai(z-h E + h+) + B+ Bi(z ~h E + h+) 

+ cos9 n \A~ Ai(z — h E + h~) + B~ Bi(z -h E + h, 

(34b) 

The coefficients a En , a En+1 , d a En , d En+1 , A+, A n , 
B+, and B~ in expressions (|32|) — (|34j) are found by impos- 
ing the continuity conditions on the wave function and its 
first derivative at the cavity interfaces. 



3 Induced emission probability 

In this section, we derive a formula for the induced emis- 
sion probability of a photon inside the cavity. To this end, 
we will use simple physical arguments rather than a rig- 
orous treatment based on wavepackets like in Ref. [29] as 
these two approaches give the same result. According to 
Eqs. (|32|) - (|33p a monoenergetic excited atom coming up- 
wards with unit amplitude upon the cavity that contains 
n photons is found to move back downwards below the 
cavity in the \a) state with amplitude d E n and similarly 
in the \b) state with amplitude d b E n+1 . As the probability 
current density of the downward wave does not depend on 
the energy (see Eq. (|21l) ). the induced emission probability 
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T , cra {n) of a photon inside the cavity containing initially 
n photons is simply given by 

,2 



and one has 



\rl b I 
\ a E,n+l | 



\ a E,n+l | 



\rl a 

\ a E,n\ 



= 1. 



(35) 



(36) 



As shown in Appendix A, this probability can be ex- 
pressed as a function of the sole wave function component 
Lp a E of the excited atom inside the cavity through the 
relation 

2 

Pem(n) = TT 2 hl ' 



u(z) <p% n (z) Ai(z - h E + 6) dz 



valid for any mode function u(z). This relation will help 
us giving a physical interpretation to our results. 

Another relation involving the induced emission prob- 
ability follows from a certain symmetry of Hamiltonian 
([1]). Denoting the induced emission probability for a total 
energy E and a detuning 5 by V em (h El S), it can be easily 
proven that for any mode function u(z) we have 

V m (h E , ~S) = V cm (h E + ~5,~S). (38) 

We study hereafter the properties of the induced 
emission probability V cul (n) depending on whether the 
atomic kinetic energy (or the lack of this energy) at 
the cavity interfaces is much higher than the interaction 
energy (\h E \ 3> h- lnt and \h E — L\ 3> h- mt ) or much lower. 
We call each case the classical and quantum regime, 
respectively. 



3.1 Classical regime : Rabi limit 

At resonance (uj = luq), it has been shown in Ref. [29j 
that the induced emission probability ()35[) reduces, in the 
classical regime {\h E \ ^> /lint and \h E — L\ 3> h int ), to the 
well known Rabi formula 1321 



where r is the classical transit time of an atom of en- 
ergy E through the cavity, indicating that way that the 
quantization of the atomic motion is unnecessary in this 
particular regime. The formula generalizing Eq. (|39| in 
the presence of a detuning 5 can be easily derived using 
the classical results of the interaction of a two-level atom 
with a quantized electromagnetic field |32j and is given by 



Pem(n) 



On 



Tl ST 

COS — %— COS 



8 • ^'„ T 2 ■ ST 

7§r sin — sm 



sm 
10 



2 O'ti 



L < h E , 

h E ^ L, 
h E < 0, 

(40) 



where n and T2 are the classical interaction times for 
< h E < L and L < h E , respectively. For h E < 0, 
the atom does not reach the cavity and no induced emis- 
sion process can occur (the classical interaction time r is 
null). For < h E < L, the atom turns back inside the 
cavity and the interaction time t\ does not depend on the 
cavity length L. For h E > L, the atom passes twice inside 
the cavity (upwards and downwards) and the interaction 
time T2 is L dependent. We will denote the time which 
passed between the first and the second interaction by T. 
In terms of dimcnsionless variables, these various times 
are given by 



hn - \ He - L 



f)' n n=2JS* + 4hl 



ST = 2S\ h 



>E 



(41) 

(42) 
(43) 



We illustrate in Figs.[3]Ja) and[3fb) the induced emis- 
sion probability Vem = Vemin) with respect to the cav- 
ity length for h E 3> /lint) respectively with and without 
quantization of the atomic motion (Eqs. (|35|) and (j40]l 
respectively). We observe that both results are in good 
agreement, except in the region where L w h E . In this 
region, \h E — L\ is not much greater than h- lnt and we are 
out of the validity domain of the approximated expression 
PO")) . The small additional oscillations observed in the in- 
duced emission probability curve where L ps h E are pure 
quantum effects originating from the quantization of the 
atomic motion. These effects arise when the cavity length 
fits nearly the classical turning points of the atoms in the 
gravitational field. In this case, the potential barrier and 
well Vn(z) play a significant role as the kinetic energy of 
the atoms is very small exactly where the potential energy 
exhibits strong variations. 



(39) 3.2 Quantum regime : mazer limit 



In the quantum regime, the induced emission probability 
■pem exhibits a completely different behavior. This is first 
illustrated in Figs. [4] and [5] that show V cm as a function 
of the total atomic energy and the cavity length, respec- 
tively, for nonvanishing detunings. Contrary to the classi- 
cal regime, narrow resonances are presently observed, even 
in the negative energy domain [E < 0, i.e., h E < 0, see 
Eq. where classically the atom does not reach the cav- 
ity and cannot therefore emit any photon in there. The 
existence of negative energy resonances reflects the fact 
that the atom can still interact with the cavity field for 
h E < because of the tunnel effect. Indeed, we have shown 
in that at resonance (6 = 0) the total potential felt 
by the atom in interaction with the cavity possesses qua- 
sibound states of negative energy. Whenever such a state 
energy matches the atomic E value, the atom can enter 
the cavity by tunnel effect and emit a photon in there. In 
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Fig. 3. Induced emission probability "P cm with respect to the cavity length (a) with and (b) without quantization of the atomic 
motion (classical regime: Ke — 1000, hint = 1 and S = —0.2). 
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Fig. 4. Induced emission probability with respect to the total Fig. 5. Induced emission probability with respect to the cavity 



atomic energy for L = 10, hint = 10 and n = 0. 



length for hE = 1, h- m t = 100 and n = 0. 



the presence of a detuning, the emission probability reso- 
nances still exist but with a modified position, amplitude, 
and width. 

From a physical point of view, the behavior of the sta- 
tionary wave functions <p E n (z) and <p% n+ i(z) [Eqs. (f32j) 
and (|33[) ] differs strongly depending on whether the value 
of the total atomic energy E gives rise or not to a reso- 
nance in the emissio n prob abi lity cu rve [see Fig. 0] . This 
is illustrated in Figs. 6(a) and |6(b) that show the atomic 
probability densities \<Pe n (z)\ an d \<Pe n +\{z)\ 2 inside 
and outside the cavity (grey and white area respectively). 
In the first case, the atom is launched upwards to the 



cavity with a total energy He = —1.90991 matching a 
resonance in the emission probability curve [see Fig. [4] 
and the stationary wave functions ip E n (z) and cp E n+1 (z) 
take a significant value inside the cavity. Consequently, a 
strong interaction between the atom and the cavity field 
can occur and the induced emission probability is signifi- 
cant, according to Eq. (|3"T)) . When considering a realistic 
atomic wave packet, this effect is of course restricted to the 
components of the wave packet for which h e falls within 
the width of the resonance. Inversely, if the to tal at omic 
energy does not match a resonance energy (Fig. 6(b) ), the 
stationary wave functions are almost zero inside the cav- 
ity and no atom-field interaction can take place. In this 
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-10 10 



(a) h E = -1.90991 [Pcm ^ 0.35] 




-10 



(b) h E = -2.2 [T em ~ 0] 

Fig. 6. Atomic probability densities \<Pe n(z)\ 2 and 
\tp b E, n +i(z)\ 2 for L = 10, hmt = 10, 5 = -1, and n = 0. The 
grey area shows the cavity region. 

case, the induced emission probability drops down to zero 
according to Eq. ([57)1 . 

3.2.1 Peak position 

Figure [JJ shows the induced emission probability for a pos- 
itive, a negative and a null detuning. We observe that for 
positive detunings, the emission probability peaks move 
towards increasing energies while for negative ones it is 
the opposite. According to Eq. ([55]) . peaks for opposite 
detunings ±5 are separated in energy by S (0.5 in the case 
of Fig. [7J. From Fig. [7J and further calculations, it turns 
out that the peak energy shift caused by a detuning ±S is 
in very good approximation given by ±6/2. 

3.2.2 Peak amplitude and width 

Figure [5] shows how the peak amplitude evolves with re- 
spect to the detuning for various fixed values of the to- 
tal atomic energy. In contrast to the classical regime [see 
Eq. (HO])] , the curves present a strong asymmetry with re- 
spect to the sign of the detuning. This results from the 



l 



0.75 



0.25 



il 



5 = -0.5 5 = S = 0.5 






1 i 
— 0.5-» 


! ! 



-3 -2 -1 



Fig. 7. Induced emission probability with respect to He for 
feint = 10, n = 0, and L > h E + feint- 

atomic energy variation hS when the atom emits a pho- 
ton inside the cavity (see Sec. I2.2[) . The smaller the total 
atomic energy is, the more pronounced this asymmetry 
is. The induced emission probability tends very rapidly to 
zero for positive detunings, in contrast to what happens 
for negative detunings. This behavior is already present in 
Eq. (|37[) where the Airy Ai function decreases much faster 
for large positive arguments than for large negative argu- 
ments as can be seen from the asymptotic expansions |30j 

Ai(-z) = ^ sin + \) + 0(z- 7 / 4 ) (44) 
Ai(z) = z~ 1/4 e-t + 0(z- 7 / 4 e-«) (45) 

valid for \z\ 3> 1 and where £ = §^ 3 ^ 2 - 

This asymmetry was also present in the absence of 
gravity [16] [horizontal geometry, see Fig.[T]. In Ref. [16], 
it had been shown that the induced emission probability 
is strictly null for h,5 > E. This results merely from en- 
ergy conservation. If the initial atomic kinetic energy E is 
lower than hS, the transition |a, n) —> \b, n+1) cannot take 
place (as it would remove HS from the kinetic energy) and 
no photon can be emitted inside the cavity. In this case 
the emission process is completely blocked. In the pres- 
ence of gravity, there doesn't exist such a threshold. For 
instance, we can see in Fig.[8]that the peak amplitude for 

= 3 takes significant values beyond hg — 3. This is not 
really surprising since quantum mechanically a particle of 
energy E moving freely in the gravity field is described by 
the wave function ip(z) = NAi(z — He) and has a non- 
vanishing probability of being found in the nonclassical 
region z > He- This probability decreases however very 
fast as a function of the height z. This enables us to un- 
derstand qualitatively why the peak amplitude decreases 
very rapidly but not sharply (no threshold as in the case 
of the horizontal mazer) as a function of the detuning for 
5 > h E - 
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Fig. 8. Amplitude A of the resonances in the induced emission 
probability with respect to 8 for h ln t = 100 and n = 0. Circles 
and crosses are the values computed from (|46|l for Ke — 3 and 
He = 1 respectively. 




-6 -3 3 6 



Fig. 10. Induced emission probability with respect to the cav- 
ity length for Ke = — 1, hint = 10 and n = 0. 

3.2.3 Cavity length effects 

In the classical regime and in the positive energy domain, 
the induced emission probability does not depend on the 
cavity length L as soon as this length is greater than He 
(see Eq. (P4TJ10 . In the quantum regime, this statement 
must be revised and holds only when L > He + This 
is illustrated in Fig. [TU] which shows the induced emis- 
sion probability "P em with respect to the cavity length for 
He = — 1 and h- m t = 10. For hs = —1, V em starts to be 
constant only beyond L — He + h+ ~ 8.512. 



3.2.4 Detuning sensitiveness 



Fig. 9. Amplitude .4 of the resonances in the induced emission 
probability with respect to Ke for hi nt = 100 and n = 0. 



Figure [H] shows the peak amplitude with respect to the 
total atomic energy for two opposite values of the detuning 
5 = ±2. According to Eq. (|38p. the curves have the same 
shape and are merely shifted in energy by He = 2. 

We have observed that the peak amplitude only 
slightly depends on the interaction height hint and that 
it is well approximated by the following simple analytical 
formula 



Figure [IT] displays the induced emission probability with 
respect to the detuning. For realistic experimental pa- 
rameters and rubidium atoms, these resonances may 
even become extremely narrow. Their width amounts only 
0.3 Hz for Ke = —2.5 and h- mt = 100. This could define 
very useful metrology devices (atomic clocks for example) 
based on a single cavity and with better performances 
than what is usually obtained in the well known Ramsey 
configuration [33] , 



3.3 Sine and Gaussian mode functions 



A 



U(-* s )D'({-ii)-D({-ij!)U'(4j) (46) 



In Fig. [51 we compare the peak amplitude derived from 
Eq. (|35)l and computed from Eq. (|46|). respectively. For 
Yie = 3 and higher values of the energy, a good agreement 
is observed. However, this agreement lowers as the energy 
decreases. 

Concerning the peak width, we have observed that 
negative detunings increase the peak width while positive 
ones decrease it (see for example Fig. [5]). 



The narrow resonances of the induced emission probabil- 
ity observed in the quantum regime are not restricted to 
the case of the mesa mode investigated so far here. Us- 
ing numerical procedures (see Appendix B) , we have com- 
puted the solutions (p E n {z) and ip b E n+1 (z) of the station- 
ary Schrodinger equation p4[) for sine and gaussian mode 
functions [u(z) — sin(7rz/L) for < z < L, elsewhere 
and u(z) — e~^~ L / 2 )~ '/ 2ry with a = L/s/2, respectively]. 
We show in Fig. [12] the related induced emission probabil- 
ity Vcm with respect to the detuning. Again, we observe 
sharp resonances whose width decreases with increasing 
interaction height (h- lnt ). 



John Martin, Thierry Bastin: Detuning effects in the vertical cold-atom micromaser 



9 




Fig. 11. Induced emission probability with respect to the de- 
tuning for Ke = —2.5, hint = 100, n = 0, and L ^> Ke + hi nt . 
For rubidium atoms, the central peak has a width of only 
0.3 Hz. 




Fig. 12. Induced emission probability with respect to the de- 
tuning for a sine mode with hint = 100, L — 10 (solid line) 
and for a gaussian mode with hint = 10, L = 1 (dashed line) 
(He = —2 and n = in both cases). 



As any mode function will give rise to such resonances 
(provided hint is large enough), the resonances in the in- 
duced emission probability will be a very general feature 
of the vertical micromaser in the cold atom regime, in 
contrast to what happens with the horizontal micromaser Appendix A 
where the induced emission probability characteristics are 
much more mode dependent [llj . This defines a very in- 
teresting property as it should help the experimenters to 
observe the particular behavior of the induced emission 
probability in the quantum regime, whatever the exact 
mode function of the cavity. Another nice illustration of 
the symmetry relation is displayed in Fig. [T31 this 
time for a sine mode function. 



Fig. 13. Induced emission probability with respect to the total 
atomic energy for a sine mode and hint = 10, L = 10 and n = 0. 



presence of a detuning have been discussed both in the 
classical and the quantum regimes. The classical results 
arc well recovered for hot atoms for which the quantization 
of the atomic motion is unnecessary. We have shown that 
atoms which classically would not reach the interaction 
region are able to emit a photon inside the cavity. In the 
quantum regime, the mazer properties are not symmetric 
with respect to the sign of the detuning and the system 
exhibits a sharp response as a function of the detuning. 
The induced emission probability has also been computed 
for sine and gaussian mode functions (see Appendix B for 
the numerical procedure we have developed). 
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In this appendix, we derive Eq. (|37[) which expresses the 
induced emission probability as a function of the station- 
ary wave function of the excited atom inside the cavity. To 
this end, let us recall the equations of motion (fH)) for the 
wave function components (f% n and (f b E n+1 (see Eq. (|T0|) 



4 Summary 

In this paper we have presented the quantum theory of the 
vertical mazer for a non resonant atom-field interaction. 
We have obtained analytical expressions for the wave func- 
tion in the special case of a constant cavity mode function. 
The properties of the induced emission probability in the 



^2 - z + h E J <p%, n (z) = h int u{z) <p Ein+1 (z), (A.la) 
d 2 

— -Z + flE-5) <P E ,n+l(z) = h intu(z) 

(A.lb) 

To solve the problem numerically, the domain is re- 
stricted to [z m i n ,5 max ]. Outside this domain, the mode 
function u(z) is regarded as null. 
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The outgoing Green function of the stationary 
Schrodinger equation of a particle of energy E in the grav- 
itational field, satisfying 



-j^ -z + h E )G E {z\z')=5{z-z'), (A.2) 



is given by [34] 

G E (z\z') = -inD(a E ^) Ai(a E ,+) (A.3) 

with 

- (z + z')±\z-z'\ 
a E ,± = -h E + - '-^ -■ (A.4) 

On the one hand, considering the right-hand side term 
of Eq. (|A.lbj) as a source term and using Green function 
(fA~3l) . the wave function component ip E n +i(^) i s given at 

Z = Zmin by 

^,„ + l(Zmin) = -™h\nt D(z min - h E + 5) X 



+ OC 



i(z) <p a EM (z) Ai(z - h E + 6) dz (A. 5) 



since by definition u(z) = for z < z m i n . On the other 
hand, ip E „ +1 (,z), which is continuous at z — z m j n , is also 
given by Eq. (fT9|) . Identification with Eq. (|A.5I) then yields 



+ OC 



d E ,n+i = -whim / u(z) (p a E ^ n {z) Ai(z - h E + 5) dz 

J —OO 

(A.6) 

whose absolute value squared gives Eq. (f3"T)) . 



Appendix B 

In this appendix, we describe the method we have de- 
veloped to solve numerically the system (|A.1[) of coupled 
differential equations with imposed boundary conditions. 



B.l Boundary conditions 

The continuity conditions of the wave function compo- 
nents (fl8)) and (|19p and their first derivatives at the cavity 
interfaces (5 m ; n and z max ) read 



VB,ra(^min) 



dip a E 



dz 



V(z mia -h E ) +d a En B(-2 



-(Zmin) = U'(z min - h E ) + d a E n D 
<A,n+l(Zrmn) = d En+1 D(z min + 5— h E ) 

^If^^-in) = D'(*mi„ +8-h E ) 



a ^"z min ~h E ) 



and 



n ( ^max ) 

dip 



a E , n Ai(i 



&l (^max) 



a% n Ai (z max - /ib) 



_ x Ai(z max + (5 - 

-1 (^max H~ 



(B.l) 



(B.2) 



B.2 Resolution scheme 

We first introduce the vector <p (z) = 

(<P%,n(x)> ( P b E,n+i(*)) • Now : let ^l(z) and ^ 2 (z) 
be two vectors satisfying (|14p with the boundary 
conditions (BC) (at z = z max ) 



(!) 



Ai(i 



dz 



Ai'(i 



^E.n+l (^roax) — 



(B.3) 



d^ 



dz 



'(^max) 



and 



Vl;,n(^max) 







-2*2- (z ) — f) 
^.n+lC^max) = Ai(i 



, dy dk" + 1 (^max) = Ai'(z max + <5 - /i_e) 



(B.4) 



The most general solution consistent with BC (|B.2[) is then 
given by 

"bright (z) = a a E n Tpi{z) + a h E n+l ~$z(z) (B.5) 

where a E n and a b E n+1 are constants. 

Next, let ips(z), ipi{z) and Tp$(z) be three vectors 
satisfying (I14| with the boundary conditions (at z — 5 m in) 



(3) 



n (-^min) 
^min j 



dz 



(4) 



V^n+lC^min) 



dz ^minj 



U(z m in - ft E ) 

U'(z min - h E ) 






(B.6) 



dz \ z mm) 



D(i 

D'( 






h E ) 



(B.7) 



and 



(5) 



' ¥>£, n (Zmin) = 
^,n+l(^mi„) = D(z n 



dz l^mmj 



h 8 - h E ) 



(B. 



The most general solution consistent with BC (|B.1[) is 
then given by 

~ftieft(z) = + d a E n + d b E n+1 ~0 5 (z) (B.9) 
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where d% n et d E n+1 are constants. 



In order to determine the value of the constants a a E n , 



d a 



and d b E n+1 appearing in 



R5I) and (IR9 



we have to solve the system of equations formed by the 
continuity conditions of the wave function and its first 
derivative at an intermediate point 5q G [z m i-n,Zr, 



^mm > '-max J 



Vleft(^o) 
^l'eft( 5 o) 



"fright (So) 
^r'ight( 5 0) 



(B.10) 



Once this system has been solved, the value of the am- 
plitude d E n , j is known and the induced emission proba- 
bility can be computed by use of Eq. 
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